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ABSTRACT 

We have obtained deep, long-slit spectroscopy along the major axis of NGC 6166, the cD galaxy in 
the cluster Abell 2199, in order to measure the kinematics of intracluster starlight at large radii. The 
velocity dispersion initially decreases from the central value of 300 kms -1 , to 200 kms -1 within a few 
kpc, and then steadily rises to 660 kms -1 at a radius of 60 kpc (Ho=75 kms _1 Mpc _1 , f2 m =0.3, f2A=0.7), 
nearly reaching the velocity dispersion of the cluster ((TA2199 = 775 ± 50 kms -1 ). These data suggest 
that the stars in the halo of the cD trace the potential of the cluster and that the kinematics of these 
intracluster stars can be used to constrain the mass profile of the cluster. In addition, we find evidence 
for systematic rotation (V/crssO.3) in the intracluster stars beyond 20 kpc. Such rotation is not seen in 
the kinematics of the cluster members. 

The surface brightness and velocity dispersion profiles can be fit using a single component mass model 
only by making unphysical assumptions about the level of anisotropy for both the stars in the cD galaxy 
and for the kinematics of the galaxies in the cluster. Two-component mass models for the cD and its halo 
are subsequently explored using the kinematics of known cluster members as an additional constraint on 
the total enclosed mass beyond the extent of the stellar kinematics. Under the assumption of isotropy, 
the observed major-axis kinematics can be reproduced only if the halo, parameterized by a generalized- 
NFW profile, has a soft core, i.e., a < 1 (a generalized- NFW halo with a = 1 is excluded due to low 
implied stellar mass-to-light ratios). This result is inconsistent with the predictions of current iV-body 
simulations for dark matter halos. 

To test the consistency of our halo profiles with those derived from strong lensing measurements in 
intermediate redshift clusters, we calculate the critical radii for tangential arcs, assuming that our best- 
fit mass models for the Abell 2199 were placed at cosmological redshifts between 0.2 < z < 0.5. The 
calculated critical radii for our best-fit two-component isotropic models range from 5" to 40", depending 
on the assumed source redshift, consistent with the radii for gravitational arcs observed in intermediate 
redshift clusters. 

We also present the results of Monte Carlo simulations testing the general reliability of velocity dis- 
persion measurements in the regime of low S/N and large intrinsic Doppler broadening. 

Subject headings: galaxies: elliptical and lenticular, cD, galaxies: kinematics and dynamics, galaxies: 
stellar content, galaxies: clusters: individual (Abell 2199), galaxies: individual (NGC 
6166), (cosmology:) dark matter 
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Fig. 1. — Mass profiles for generalized NFW dark matter halos with a £ 0, 1, 1.5 (solid, dashed, dot-dashed, respectively). These mass profiles 
have been renormalized to have identical mass at 1 Mpc. The left and right panels are identical, except for the linear and logarithmic stretch 
in radius. Note the range of projected radii, R < 100 kpc, over which one would observe significant differences. 



Spergel & Steinhardt 2000; Firmani et al. 2001; Colin 
et al. 2000). Halos of self-interacting dark matter have 
soft, uniform density cores (a < 1; Yoshida et al. 2000; 
Dave et al. 2001), but their low cllipticities are considered 
problematic (e.g., Miralda-Escude 2002). 

Some CDM simulations indicate that a itself may be a 
function of halo mass, with cluster-sized halos tending to 
a = 1 and galaxy-sized halos tending to a = 1.5 (Jing & 
Suto 2000, though sec Klypin et al. (2001)). Such re- 
sults suggest that the density profiles of dark matter halos 
are quite sensitive to their histories of formation, accre- 
tion, and merging, as well as to the initial conditions of 
collapse (Nusser & Sheth 1999), and these notions are 
tested explicitly by Weinberg (2001), who finds that noisy 
evolutionary processes can drive halo profiles to a uniform 
shape with soft cores. This result appears similar to the 
profile derived by Iliev & Shapiro (2001), wherein the col- 
lapse of spherical top-hat density perturbations was shown 
to produce truncated isothermal spheres, i.e., a — 0. The 
solutions of Iliev & Shapiro (2001) and Weinberg (2001) 
may naturally produce soft cores because of the bound- 
ary condition of finite phase-space densities at r = 0. This 
restriction on the solutions to the Lame-Emden or Fokker- 
Planck equations (Iliev & Shapiro 2001; Weinberg 2001) 
runs counter to the basic assumption in fitting Eq. 1 to 
TV-body profiles, wherein the central phase-space densities 
are allowed to diverge (Taylor & Navarro 2001). 

While such analytic work is providing insight into the 
physics of dark matter halo evolution, there are few obser- 
vational constraints on cluster scales to test the 1 < a < 
1.5 models routinely assembled by the simulations. For the 
massive dark matter halos of rich clusters, this range in a 
yields models that are quite similar outside of < 100 kpc, 
where most cluster- related observations have weight. Fig- 
ure 1 shows the mass profiles of example generalized- NFW 
mass profiles for a e {0,1, 1.5}. The two panels show the 
mass profiles using linear and logarithmic stretches in clus- 
ter radius and illustrate the need for observations within 
the inner < 100 kpc in order to measure the inner slope of 
the dark matter halo profile (see also Klypin et al. 2001). 

Observational tests at these small radii are scarce. The 
determination of the mass profile from kinematics of clus- 



ter galaxies within this region has been difficult because 
there are generally few members, van der Marel et al. 
(2000) attempted to solve this problem by combining the 
rcdshifts from the CNOC1 cluster survey (Yee, Ellingson, 
& Carlberg 1996) into a single fiducial cluster. They 
argued that the best-fit mass profiles were described by 
awl, but their analysis did not take into account the con- 
tribution of the stellar mass of any brightest cluster galax- 
ies (BCG) to the total mass profile of the fiducial cluster. 
An accurate determination of mass profiles of the inner 100 
kpc for clusters from X-ray observations has been ham- 
pered to date by poor spatial resolution, and by the pres- 
ence of cooling flows and shocked regions. Even as X-ray 
observations improve, uncertainties about the mass con- 
tributions from BCGs will persist (see Arabadjis, Bautz, 
& Garmire 2001). Still another method, that of using the 
kinematics of planetary nebulae and/or globular clusters 
(e.g., Cohen & Ryzhov 1997; Cote et al. 2001), requires 
painstaking identification and spectroscopy of 20-200 dis- 
crete objects, and Hcrnquist & Boltc (1993) show that the 
resulting inferred mass profiles are subject to large (> fac- 
tors of 2) uncertainties with such small samples. Strong 
gravitational lensing measurements can allow one to in- 
fer the total surface mass density within the Einstein ra- 
dius but these constraints are susceptible to degeneracies 
(Gorenstein et al. 1988) and to uncertainties in the details 
of the potential (Bartelmann et al. 1995). The mass pro- 
file can be better-constrained when there are several arcs, 
but the brightest cluster galaxies may themselves domi- 
nate the shapes of the total surface mass density profile, 
as may be the case in MS2137-2353 (Hammer et al. 1997). 
Weak lensing measurements are even less effective at con- 
straining the inner mass profile, given the typical number 
densities of background galaxies as well as the compara- 
tively large surface areas covered by the galaxies in cluster 
cores (Hockstra et al. 1998). 

A more promising method for constraining the mass pro- 
files in the cores of rich clusters is to use the kinematics 
of stars in the extended halos of central cD galaxies ( "in- 
tracluster starlight"). Past efforts to trace the stellar ve- 
locity dispersion profile of IC 1101, in the cluster A2029, 
to tens of kpc showed that a(r) rises from 400 kms -1 to 
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Fig. 2. — A log-stretch of a B-band image of NGC 6166 with the position of the LRIS long-slit shown for the primary integrations (totaling 
7200 s). The radial bins, from which the galaxy spectra were extracted, are shown. Note the increasingly large apertures from which the 
two-dimensional spectra were summed to achieve adequate signal-to-noise ratio. A second long-slit position, exposed for 1400 s, was also 
obtained with the slit centered on the galaxy nucleus. 



600 kms 1 towards larger radii (Dressier 1979; Carter 
etal. 1985; Sembach & Tonry 1996). Of the many BCGs 
investigated over radii of 10" — 20", only IC1101 (Fisher, 
Illingworth, & Franx 1995; Kronawitter et al. 2000) and 
NGC 6166 (Carter, Bridges, & Hau 1999) have shown such 
a rise. Such results suggest that the stars in a cD's halo 
respond to the potential of the cluster. 

In this paper we report on a pilot study in which we 
targeted NGC 6166, the giant elliptical in the center of 
the rich cluster Abell 2199 (z = 0.03), with the Keck II 
10m telescope and the Low Resolution Imaging Spectro- 
graph (Oke et al. 1995, ; LRIS). Several factors make this 
galaxy an ideal choice: (1) it is at rest with respect to the 
center of the cluster potential (Zabludoff et al. 1990); (2) 
the cluster itself appears relaxed in X-rays; (3) the pres- 
ence of a cooling flow in the cluster core suggests little 
or no recent merger activity in the cluster center; (4) the 
useful features in its spectrum (e.g., the G-band) do not 
overlap any strong night-sky emission lines or absorption 
lines; (5) the rcdshift allows the 170" LRIS slit to sub- 
tend 100 kpc (H = 75 kms" 1 Mpc~\ Q M = 0.3, and 
Oa = 0.7; used throughout this paper); and (6) its dec- 
lination of 39° makes it observable for most of the night 
from Mauna Kea. The telescope's large aperture makes it 
possible to obtain spectra with sufficiently high S/N ratios 
in the low surface brightness halo of the galaxy while using 
a narrow slit-width (both to achieve moderate resolution 
and to exclude interloping galaxies). 

We structure the paper as follows. In §2 we discuss the 
observations and data reduction. The stellar mass density 
profile is constrained in §3 through the use of the B- and R- 



band surface brightness profiles, and the major-axis kine- 
matics of NGC 6166 are derived in §4. In §5 we discuss 
the reliability of velocity dispersion measurements in the 
regime of low S/N and large intrinsic Doppler broadening 
using Monte Carlo simulations. Section 6 discusses the 
fit of simple mass models to the kinematics of NGC 6166 
and Abell 2199, as constrained by the kinematics of clus- 
ter members. The implications of our best-fit mass models 
are discussed in §7 and our conclusions are summarized in 
§8- 

2. OBSERVATIONS AND DATA REDUCTION 

The primary goal of this work is the accurate measure- 
ment of the second moment (a) of the line-of-sight velocity 
distribution (LOSVD) as a function of position along the 
major-axis of NGC 6166. By combining the a profile with 
broad-band imaging, we can obtain constraints on the to- 
tal mass-to-light {M/L) profile of the galaxy and the in- 
ner portions of the surrounding cluster. In this section we 
describe the B- and i?-band imaging of NGC 6166 and 
its halo, the long-slit spectroscopy, and the processing re- 
quired to derive the chief products: the surface brightness 
profile, n(R), and the projected line-of-sight velocity dis- 
persion profile, o~(R). Using the cosmological parameters 
given in §1, the scale on the sky is 0.56 kpc/arcsec. 

2.1. The Photometry 

Three £>-band images of NGC 6166, totaling 1870 sec- 
onds, and six i?-band images, totaling 1280 seconds, were 
obtained with LRIS during the nights of 30-31 May 1998. 
The images span 6x8 square arcmin, corresponding to 
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Fig. 3. — (left) The mean R- (solid) and B-band (dashed) surface brightness profile of NGC 6166, derived from isophote fitting, as a 
function of the major-axis of the isophote. (middle) The R- and B-band surface brightness profile of NGC 6166, derived along the position of 
the LRIS long-slit. The dotted line in both panels indicates the King model that best fits the B-band profile as measured along the position 
of the long-slit (see §6. The fit of the King profile was restricted to radii 2"5 < r < 125" (1.2 kpc < r < 70 kpc) and satisfactorily matches 
the profiles determined from the isophotes and from the long-slit position, (right) The B — R color gradients derived from the isophote fitting 
(solid line) and the long-slit position (dashed line). Note that the cD becomes bluer with projected radius. The color gradient is indicative 
of a radial variation in the metallicities and/or ages of the stellar populations and implies that the stellar M/L ratio is also varying with 
projected radius. 



a distance along the major axis of 170 kpc. The seeing 
during the run was approximately 0.8 arcsec (FWHM). A 
subsection of the i?-band image is shown in Figure 2. 

2.1.1. The Surface Brightness Profile 

The sky background was determined for each image us- 
ing the corners of the image furthest from the galaxy, and 
the sky brightness was then subtracted from the entire im- 
age (in both bands). While errors in the sky will severely 
affect the derived surface brightness profile of the galaxy 
at radii approaching the extent of the imaging field-of- 
view, the accuracy in our estimate of the sky is suitable 
for deriving the surface brightness profile over the extent 
of the LRIS long-slit (~3'). The B- and .R-band sky bright- 
nesses were approximately HB,sky = 21.9 mag/arcsec 2 and 
HR,sk y = 20.6 mag/arcsec 2 , respectively. 

The fitting of elliptical isophotes was performed using 
the ellipse package of STSDAS in IRAF. The detection 
of 15 interloping galaxies necessitated modeling (using the 
ellipse package) and subtraction of their light from the 
images in both bands. Isophotes were fit to these galax- 
ies, as well as to the cD, and iteratively subtracted from 
the image. The envelope of one interloper intersects the 
major-axis between 2 kpc and 14 kpc at a mean level of 
about 15% of the cD surface brightness. The residuals 
from the subtraction of this particular interloper amounted 
to -5% of the cD surface brightness at those radii along 
the major-axis. Such residuals were efficiently rejected by 
the sampling criteria in the final ellipse-fitting of the cD. 

2.1.2. The Calibration 

We calibrated the -B-band photometry using secondary 
standard stars in the Galactic globular cluster M92, with 
an error in the photometric zero-point of ±0.016 mag. Cal- 
ibrating the i?-band photometry was problematic because 
the LRIS filter wheel had not properly inserted the (John- 
son) R filter into the optical path for the M92 calibra- 
tion image. As a result we used the Cousins i?-band sur- 
face brightness profile from Postman & Lauer (1995) to 
calibrate our data. The Postman & Lauer (1995) profile 



was first transformed from Cousins R to Johnson R using 
Rj — Rq = —0.11 (a good approximation for the spectral 
energy distributions of old E/SOs; (Fukugita, Shimasaku, 
& Ichikawa 1995). We then determined the photometric 
offset between the LRIS profile and the transformed Post- 
man & Lauer (1995) profile over the radial range of 5 to 
40 arcsec. The rms of this offset was 0.017 mag. 

The i?-band extinction to NGC 6166, as estimated by 
Schlegel et al. (1998), is A B = 0.05 mag (A B = 0.00 mag 
from Burstein & Heiles 1982). Because this prediction is 
quite low and irrelevant for the determination of the mass 
profile we have made no correction for the foreground ex- 
tinction (it only affects the deduced stellar Mj L ratios at a 
level of < 5%). The final B- and i?-band surface brightness 
profiles are shown in the left panel of Figure 3. 

2.1.3. The Surface Brightness Profile along the LRIS 

Long-Slit 

While isophote fitting provides the maximum signal-to- 
noisc ratio in the surface brightness profile at a given ra- 
dius, our kinematic profile will be derived along a fixed 
position angle (the major-axis). Thus it is more appropri- 
ate to use the surface brightness profile along this position 
angle in modeling the mass profile and constraining the 
mass model of the cD and its halo. 

We calculated the B- and i?-band slit profiles from eight 
columns of the CCD image centered on the galaxy nu- 
cleus. The long-slit intensities were derived by averaging 
the eight columns of the slit position in order to produce 
a single intensity per pixel for each row along the slit. 
These B- and i?-band long-slit surface brightness profiles 
are shown in the middle panel of Figure 3. The profiles in 
the right- and left-hand panels are virtually identical be- 
cause there was very little twisting of the isophotes with 
radius. Thus our mass modeling will be insensitive to the 
choice of the long-slit or isophote profile. For the remain- 
der of the paper, we use the profile derived at the position 
of the LRIS long-slit. 

2.1.4. The Color Gradient m NGC 6166 
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Fig. 4. — The instrumental resolution for the LRIS configuration used for the work presented in this paper. The resolution was determined 
from the widths of the emission lines from the mercury, krypton, and neon lamps. 



The right-hand panel of Figure 3 shows the B — R color 
gradients derived from both the isophote fitting (solid line) 
and for the long-slit position (dashed line) . The color tran- 
sitions smoothly from B—R = 1.6 mag in the inner regions 
of the cD to B-R= 1.4 mag at r = 20 kpc to B-R = 1.2 
mag at r = 50 kpc. These are colors which are typical for 
E/SO, Sa, and Sbc galaxies, respectively (Fukugita, Shi- 
masaku, & Ichikawa 1995). We find that the gradient is 
well fit by the second-order polynomial: 

(B-R) = 1.58-0.10[logr(kpc)] 2 (2) 
with an rms scatter of 0.018 mag. 

In §3 we derive structural parameters from the surface 
photometry and discuss the use of the color gradient to 
estimate possible variations in stellar M/L with projected 
radius. 

2.2. The Spectroscopy 

The primary spectroscopic data consist of four 40 
minute exposures of the galaxy, flanked by five 40 minute 
offset-sky exposures (taken before and after the on-galaxy 
exposures). We used the 1"5 wide long-slit with the 600 
grooves mm -1 grating, blazed at 5000 A. The spectral cov- 
erage extends from 3800A to 6400A. The resolution of ~7A 
FWHM is equivalent to CT; ns t. ~ 200 km/s, as shown in Fig- 
ure 4, more than adequate to resolve the velocity disper- 
sions in the cD and its halo. While the spectral coverage 
is suitable for deriving absorption line velocity dispersions 
from several features, e.g., the Mg b, G-band, or Ca H 
& K, we concentrated our efforts around the G-band as a 
compromise between the desires to minimize contamina- 
tion due to the sky background and to achieve adequate 
signal-to-noise ratios. 

To maximize the extent of the halo observed and still 
preserve radial profile information for NGC 6166, we 
placed one end of the long-slit at the galaxy center and 
the rest of the 170" long-slit north-east along the major 
axis at a position angle of 33 degrees (see Figure 2. Our 
imaging of NGC 6166 on the same night confirms this an- 
gle as that of the major axis to within a degree. A second 



position, at the same position angle but centered on the 
galaxy nucleus, was also obtained in 1400 sec of integration 
time. 

Because NGC 6166 fills the entire length of a long-slit 
placed long the major-axis, we obtained alternating sky 
pointings offset from the galaxy, similar to the approach 
used by Dressier (1979), to facilitate accurate removal of 
the sky background. In principle, obtaining the sky spec- 
trum at the same position on the detector as the galaxy 
spectra should minimize the effects of imperfect flat field- 
ing, errors in rectification, and uncertainties in interpolat- 
ing the sky background from the slit ends. While Sem- 
bach & Tonry (1996) used charge shuffling and chopping 
of the telescope on 5 min timescales in order to ensure 
that the offset-sky spectra were obtained at precisely the 
same physical pixel locations as the galaxy spectra, LRIS 
cannot be used exactly in this mode. Therefore we chose 
to nod between galaxy and sky on 40 min intervals, having 
to read out the detector at each position/exposure. Our 
exposure times were chosen as a compromise between (1) 
requiring the data in the low-surface brightness halo of the 
cD to be limited by the noise in the sky background and 
not by the read noise in the electronics, and (2) avoid- 
ing the time-variability of the background. Despite the 
variability of the sky spectrum, and the flexure character- 
istics of the instrument, we were able to obtain quality 
halo spectra where the surface brightness of the galaxy is 
approximately 6% of sky (see below). 

2.2.1. Rectification and Wavelength Calibration 

We used the software package Expector (Kelson et al. 
2000) to rectify and wavelength-calibrate the LRIS spec- 
troscopic data. A combination of Hg, Kr, Ne, and Ar 
lamps allowed us to rectify the spectra and to derive the 
high order terms for the wavelength calibration. The rms 
scatter about the dispersion solution of the lamp spectra 
was 0.04A (3 kms" 1 at the G band). The zero-points of 
the dispersion solutions of the galaxy spectra were refined 
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Fig. 5. — A spectrum of the Mauna Kea night sky in the blue. Note the tremendous amount of structure blueward of O I 5577A, including 
the (strongly variable) nitrogen line at 5199A, which is superposed on OH Q & R lines at ~5200A. Because most of the features are broad 
with very sharp edges, poor sky subtraction of low-surface brightness sources can artificially alter the broad absorption features intrinsic to 
a galaxy spectrum. 



using several high S/N sky emission lines (e.g., 5224A, 
5238 A, 5577A, 6300A), resulting in a final rms scatter 
of 0.6A (±40 kins" 1 at the G-band). This large resid- 
ual scatter about the sky lines is typical for LRIS data 
largely because the lines in lamp spectra are asymmetric 
compared with the sky emission lines. The asymmetries 
in the line profiles systematically vary with position across 
the CCD causing small systematic errors in the measure- 
ments of the line centers, and thus introduce a bias in the 
wavelength scale (Kelson et al. 2000). These large resid- 
uals appear as a moderate-order deviation from the lamp 
line dispersion solution and may be reduced with further 
work if one intends to utilize a very long spectral baseline 
(see Kelson et al. 2000). If the direct fitting of broad- 
ened templates to the galaxy spectrum is performed over 
a relatively short wavelength range, then a moderate rms 
in the dispersion solution manifests itself primarily as a 
zero-point error, and secondarily as a much smaller effec- 
tive rms scatter in a local dispersion solution. For our 
data and fitting range, the effective scatter is smaller by a 
factor of four (to ±10 kins" 1 ), and thus has a negligibly 
small impact on the measured velocity dispersions. 

2.2.2. Background Subtraction 

By obtaining full 2D sky spectra at offset positions from 
the cD, before and after each galaxy exposure, one should, 
in principle, be able to reconstruct the sky background at 
the time of the galaxy exposures. We modeled the tem- 
poral variability in the sky by fitting low-order polynomi- 
als to the sky values at each pixel, as functions of time 
and/or airmass. This method should have produced the 
expected 2D spectra of the sky background at the time 
and/or airmass of the galaxy observations, but such simple 
algorithms and interpolations did not allow for sufficiently 
accurate subtraction of the variable sky spectrum (even in 
the blue; see Figure 5. 

We improved the background subtraction by fitting each 
CCD row in the on-galaxy 2D exposures using Equation 8 
of Kelson et al. (2000), originally suggested for measuring 



velocity dispersions from data from which the background 
has not been subtracted. By doing so, each CCD row in 
the data was fit as a sum of (1) the high signal-to- noise sky 
spectrum (from the offset sky-exposures in the same CCD 
row), (2) a broadened velocity dispersion template, and 
(3) the low-order continuum functions. Thus, the data 
were represented by a sum of the galaxy's contribution 
with the spectrum of the night sky, scaled by and offset by 
low-amplitude and low-order polynomials. This polyno- 
mial, equivalent to < 1% errors in the sky determination, 
effectively removes the time-variability of the sky back- 
ground and any instabilities in the instrument. Additional 
broadening of the sky spectrum did not improve \ 2 so we 
infer that focus variabilities were negligible. Construction 
of the background in this way provides a high signal-to- 
noise map of the contribution of the sky spectrum to the 
data within each row. At large galactic radii, where the 
galaxy contributes very little flux, this construction of the 
row-to-row map of the background produces a noisy row- 
to-row map of galaxy flux. Therefore, the 2D map of the 
sky contribution was fit by a low order polynomial in each 
wavelength bin and this spatially smoothed map of the 
sky background was used to perform the background sub- 
traction. Based on the analysis in §5, in which we discuss 
the robustness of the a measurements, our results are not 
sensitive to any remaining uncertainties in the background 
subtraction. 

2.2.3. Extraction 

Spectra were summed along the spatial direction of the 
CCD in order to produce ID spectra with adequate signal- 
to-noise ratios (> 20 per A) for determination of velocity 
dispersions. Figure 2 illustrates the positions of the end- 
points of the spatial bins using short solid lines perpen- 
dicular to the length of the slit. In Figure 6 we show the 
total sky±galaxy spectra from the deep spectroscopic in- 
tegrations of the cD halo in four representative bins along 
the major-axis. The corresponding sky-subtracted galaxy 
spectra are shown in the right-hand panels of Figure 6. 
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Fig. 6. — The left-hand panels show the total sky + galaxy spectra in four representative bins in projected radius along the major axis of 
NGC 6166. The right-hand panels show the galaxy spectra extracted from the same radial bins. The velocity dispersion fits, shown below the 
galaxy spectra, illustrate the quality of the fitting method. The residuals from each fit are also shown and these indicate that no significant 
systematic effects remain unaccounted. 



The data have not been flux-calibrated. 

3. CONSTRAINTS ON THE STELLAR LUMINOSITY DENSITY 

In constructing mass models for NGC 6166, the density 
profile of the stellar mass component is constrained by the 
observed surface brightness profile. The surface brightness 
profile only constrains the shape of the stellar luminosity 
density profile; the stellar M/L ratio (M/L*) can only be 
constrained by the kinematics. Traditionally, photometry 
from a single bandpass is used, with the assumption of a 
spatially invariant M/L*. However, galaxies have gradi- 
ents in their stellar populations (e.g., Peletier et al. 1990; 



Franx, Illingworth, & Heckman 1989; Davies, Sadler, & 
Peletier 1993; Gonzalez 1993; Fisher, Franx, & Illing- 
worth 1996; Kuntschner 1998; Trager et al. 2000) and by 
inference must have gradients in M/L*. Fortunately, we 
have surface photometry in both B and R, and in this sec- 
tion we combine the B — R color gradient from §2.1.4 with 
simple stellar population synthesis models to estimate the 
gradients in the S-band M/L* ratio (M/Lb) and i?-band 
M/L* ratio (M/L R ). 

While the nature of the color gradient (as a spatial vari- 
ation in stellar ages and/or in stellar chemical abundances) 
remains uncertain, we can make some simple assumptions 
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to test to what extent the gradient in the stellar mass- 
to-light ratio will affect our conclusions about the profile 
of gravitational mass in the system. We explore three ex- 
treme cases: (1) no gradient in the properties of the stellar 
population; (2) a gradient in [Fe/H]; and (3) a gradient in 
stellar ages. Using the Vazdekis ct al. (1996) models, we 
have derived the correlation of (B — R) color with Mj L b 
and M/Lr at fixed age or at fixed metallicity: 
Alog(M/L B )| fixed age = 0.96A(B- R) 
A log(M/L R )\ &xed agc = 0.56A(B - R) (3) 

and 

Alog(M/L B )| fixed [Fe/H] = 2.00A(B - R) 

A log(M/L JJ )| fixed [Fe/H] = 1.60A(B - R). (4) 

Combining the above relations with the color gradient in 
NGC 6166, we find that the adoption of a constant M/Lb 
would lead one to overestimate the stellar mass densities 
at large radii by factors of 2-6. By incorporating the gradi- 
ent in stellar populations, we hope to derive (presumably) 
more accurate representations of the stellar mass density 
profile (though the results would still be dependent on the 
assumption of the age- or metallicity-dcpendcnce of the 
color gradient). 

In Table 1 we list the parameters of the best-fit de Vau- 
couleurs, King (1966), and power-law profiles, under the 
three extreme assumptions about the gradient in the stel- 
lar populations (the power-law profile is assumed to be of 
the form: p oc [1 + (r/r c ) 2 ] 7 ). All three parameterizations 
fit the data well, with rms residuals between 0.03 and 0.07 
mag. 

When fitting the model profiles to the uncorrected B- 
and i?-band surface brightness profiles, one obtains smaller 
effective and core radii in the redder bandpass, consistent 
with the stellar populations becoming bluer with increas- 
ing radius. When the B- and i?-band surface brightness 
profiles are corrected, the resulting structural parameters 
from the two bandpasses are essentially identical. Because 



we have divided our observed surface brightness profiles by 
the gradients in M/Lb and M/Lr, the corrected surface 
density profiles should more accurately reflect the stellar 
mass density profile (modulo the normalization by a con- 
stant) . 

As an aside, there are large differences in the values of 
r e obtained from the uncorrected, age- variation-corrected, 
and metallicity-variation-corrected profiles. In the case 
where the stellar population gradient of NGC 6166 (and 
presumably of other giant ellipticals) is due solely to gra- 
dients in the chemical abundances, the deduced half-mass 
radius is smaller than the B- and i?-band half-light radii 
by 45% and 30%, respectively. Such large systematic er- 
rors in the half-mass radii of elliptical galaxies may have a 
grave impact on early-type galaxy scaling relations, such 
as the fundamental plane (Faber et al. 1987; Djorgovski 
& Davis 1987), through the strongly correlated errors in 
r e and ji e and the modest scatter in the color gradients 
of ellipticals and SOs (e.g., Peletier et al. 1990). For the 
work presented here, we do not utilize the half-light or 
half-mass radii, so we offer this short discussion simply as 
a cautionary note. 

Now that we are free from the assumption of a constant 
stellar mass-to-light ratio, we can explore to what extent 
the stellar mass may or may not be sufficient to support 
the motions of the stars through the cD halo. 

4. THE MAJOR-AXIS KINEMATICS OF NGC 6166 

Spectra for the galaxy were extracted from the radial 
bins defined in Figure 2. Four represenative spectra are 
shown in the right-hand panels of Figure 6. In each of the 
radial bins, S/N > 20 per A, with the S/N reaching 120 
per A inside 15 kpc. With such S/N ratios, the data are 
quite suitable for measuring the internal motions of the 
stars in NGC 6166 and its halo out to 60 kpc. To simplify 
the derivation of the internal kinematics, we parameter- 
ize the LOSVD as a Gaussian. In the inner regions the 
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Fig. 7. — (left) The observed projected velocity dispersion profile of NGC 6166 derived from the deep, off-set long-slit pointing (filled 
circles), and from the shallow, central long-slit exposure (open circles and squares). The squares indicate the profile from the opposite side 
of the galaxy from the filled circles, with the projected radii flipped for comparison. The open and closed circles between 4 kpc < R < 12 
kpc were obtained by simultaneously fitting both the primary (NGC 6166) and secondary (interloper at +1300 km/s) component. Note that 
the velocity dispersion profile has a large positive gradient in a on both sides of the galaxy after a decline from the central value of <ro = 300 
km/s. (center) The reduced \ 2 as a function of position. Our methods provide excellent fits to the NGC 6166 spectra except in the outermost 
bin, where the galaxy surface brightness is 6% of the sky in the B-band. Between 4 kpc < r < 12 kpc the low reduced x 2 values indicate 
that the interloping galaxy at Sv = +1300 km/s was successfully removed in the fitting procedure and had negligible impact on the velocity 
dispersions measured for NGC 6166. (right) The rotation curve from the same data. Note the apparent increase in rotation in the outskirts 
of the cD. 
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Fig. 8. — (left) The level of rotational support as a function of projected radius along the major-axis (symbols as in Fig. 7. Most of the stars 
in NGC 6166 are clearly supported by random motions with an apparent increase in rotational support outside of R = 20 kpc. (center) The 
cllipticity of the B-band isophotes as a function of projected radius along the major-axis. The isophotes of NGC 6166 become increasingly 
elliptical at large radii, consistent with an increase in rotational flattening, (right) The level of rotational support as a function isophotc 
cllipticity in the S-band. The curves show the expected correlations between \V/a\ and e for rotationally supported oblate spheroids with 
(5 = 1 — cfzz/cj xx e {0,0.1,0.2,0.3,0.4,0.5} observed edge-on. Using these data we conclude that within a projected radius R < 20 kpc NGC 
6166 is not a rotationally flattened oblate spheroid. Beyond projected radii R = 20 kpc, the isophotes become more elliptical and the cD 
halo appears to become mildly anisotropic and rotationally flattened. This level of rotational flattening only amounts to < 5% of the kinetic 
energy in the halo and has little effect on our results. 



high S/N ratios would normally have allowed us to mea- 
sure the higher-order moments h 3 or /14 of the LOSVD 
(van der Marel & Franx 1993). Unfortunately (see be- 
low), an interloping object with a velocity shift of 1300 
km/s from NGC 6166 contaminates its spectrum between 
4 kpc < r < 12 kpc. As a result the LOSVDs for both 
objects were modeled by pure Gaussians, with both com- 
ponents being fit simultaneously (e.g., Kelson et al. 2000). 
At large radii in the cD halo, where accurate determina- 
tions of h 3 or hi would also be useful for constraints on 
the anisotropy, the S/N ratios are not sufficient. How- 
ever in other giant ellipticals the departures of LOSVDs 
from Gaussian are typically small, even at large radii (e.g., 
Kronawitter et al. 2000). 

For each extracted spectrum, we performed a least- 
squares fit of broadened template spectra, using stellar 
template spectra from the literature (Kelson et al. 2000) 
and a more narrow range of stellar spectra obtained dur- 
ing the observing run. Using the procedures outlined in 
Kelson et al. (2000), we use night-sky air-glow emission 
lines and calibration lamp lines to measure the instrumen- 
tal broadening. Once the instrumental broadening had 
been characterized as a function of wavelength and posi- 
tion along the CCD, the template spectra from the litera- 
ture were broadened to match the instrumental resolution 
of the galaxy spectra. While we relied primarily on the 
templates from the literature to derive the internal kine- 
matics, the template spectra obtained during the run were 
useful to verify the robustness of our measurements. The 
velocity dispersions derived using the templates observed 
during the run agreed with velocity dispersions derived us- 
ing the templates from the literature at the level of 0.5%. 

The quality of the fits to the galaxy spectra can be seen 
in the right-hand panels of Figure 6. While the inter- 
loper is clearly visible in two of the examples in the fig- 
ure, the contamination was easily removed by the fitting 
procedures and the residuals of the fits show no remain- 
ing systematic effects. The other two radial bins shown 



in the figure show the quality of the velocity dispersion 
fits at large projected radii (17.2 < R < 20.1 kpc and 
37.3 < R < 41.3 kpc). 

In the left panel of Figure 7, we use filled circles to show 
the velocity dispersion profile derived from the deep, pri- 
mary long-slit observations on one side of the galaxy (see 
Figure 2). The open points (circles and squares) show the 
data from the shallow central pointing; the open circles 
are from the same side of the galaxy as the deep pointing, 
and the open squares were derived from the opposite side 
of the galaxy. The agreement of a(R) from both sides of 
the galaxy is excellent and indicates that the interloping 
galaxy was successfully removed during the fitting of the 
deep spectroscopy. The reduced % 2 values, shown in the 
central panel, show the quality of the velocity dispersion 
fits. 

In the deep spectroscopy of the cD halo there is evidence 
for a significant level of rotation beyond R > 30 kpc, con- 
sistent with the results of Carter, Bridges, & Hau (1999) 
to R <~ 20 kpc. The right panel of Figure 7 shows the stel- 
lar rotation curve and the fraction of rotational support 
(|U/er|) as a function of projected radius is shown in the 
left-hand panel of Figure 8. 

The error bars shown in the figures represent the formal 
uncertainties only. Despite the 18% formal uncertainty 
and elevated x 2 m the last radial bin, the gradient in a is 
clearly significant and is also consistent with the profile of 
Carter, Bridges, & Hau (1999) to R < 20 kpc. Additional 
errors, such as those arising from template mismatch, are 
estimated at ~3-5% by comparing the velocity dispersions 
using templates of varying spectral types. Errors in a due 
to background subtraction are included in the estimates 
of the formal uncertainties (verified by using the y/ 2 mini- 
mization described in §5.2.4 of Kelson et al. 2000). The 
implications of the major-axis kinematics for the mass pro- 
file of NGC 6166 will be discussed in §6. 

Table 2 lists the major-axis kinematics for NGC 6166. 
The columns are as follows: (1) the central radius of the 
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Fig. 9. — Monte Carlo simulations with more than 3 X 10 4 realizations of our NGC 6166 spectra for which (1) there is no mismatch between 
the template and the artificial galaxy spectra; and (2) surface brightness 7j n > 0.05/ s k y - The upper panels show the true fractional errors 
in a as a function of "observed" velocity dispersion, c (, s . The "observed" parameters <7 (, s and the formal errors, <5 j, s , were derived by 
applying the direct fitting method of Kelson et al. (2000) to the simulated spectra. The velocity dispersions input into the simulated spectra 
are denoted by Uj n . Any systematic offsets between the points in the upper panels from zero would indicate a decrease in our ability to 
accurately measure velocity dispersions within a given range of cr obs . The dashed lines in the upper panels show the systematic error in cr obs 
as a function of S/N ratio. These diagrams indicate that our method can reliably recover velocity dispersions up to a b s < 700 km/s. For 
larger values of a, there appear to be systematic errors > + 5% for S/N < 30 per A. The lower panels show the observed formal errors as a 
function of S/N ratio in the given velocity dispersion ranges. The mean formal error as a function of S/N is shown by the solid line while the 
"true" formal errors, as derived from the scatter in the upper diagrams, is shown by a dashed lines in each of the lower panels. The agreement 
between the two curves illustrates that our formal errors arc reliable estimates of the random error in the measurement. 



radial bin, in kpc; (2) the length of the radial bin, in kpc, 
(3) a, the velocity dispersion, in kms -1 ; (4) 8 a , the for- 
mal error in a, in kms -1 ; (5) V, the mean line-of-sight 
radial velocity, in kms^ 1 ; (6) <5y, the formal error in V, in 
kms -1 ; and (7) the reduced \ 2 of the fit. 

5. MONTE CARLO SIMULATIONS OF THE ACCURACY OF 
VELOCITY DISPERSION MEASUREMENTS 

Accurate measurement of velocity dispersions, i.e., ab- 
sorption line widths, from galaxy spectra requires: (1) that 
the instrumental broadening of the galaxy spectrum must 
be (approximately) less than the intrinsic Doppler broad- 
ening; and (2) that there must be statistically significant 
structure on scales comparable to the Doppler broadening. 
If either of these requirements is not satisfied, systematic 
errors may become important. The great difficulty in mea- 
suring the velocity dispersion profiles of galaxies to large 
radii, and for cD galaxies in particular, is that this last 
requirement may not be satisfied at large projected radii. 
In the cD halo, the potential may become dominated by 
dark matter and the stars may have similar kinematics 
to the cluster galaxies. Furthermore, the surface bright- 
ness of the stellar component of the cD halo is typically 
quite low (several magnitudes below the dark sky in the 
_B-band) . At such low surface brightnesses even a small er- 
ror in the background subtraction can introduce unwanted 



structure/power into the extracted galaxy spectrum. Be- 
cause the representative sky spectrum in Figure 5 has sig- 
nificant power on the scales relevant for the measurement 
of absorption line velocity widths, there exist the possi- 
bilities of over-estimated S/N ratios and of excess struc- 
ture in the spectra at a sufficient level to bias the deduced 
LOSVD. In this section we analyze simulations of our data 
to determine under what conditions such systematic effects 
become important and at what level our measurements of 
the kinematics in NGC 6166 may no longer be trustworthy. 

Because our data span a wide range of mea- 
sured velocity dispersions over a broad range of 
surface brightnesses, we constructed 4680 artificial 
galaxy spectra with intrinsic broadening of <jj n € 
{100, 200, 300, 400, 500, 600, 700, 800, 900} km/s. We used 
the 10 Gyr-old, solar metallicity, blue SED from Vazdckis 
(1999) as the underlying spectrum of the artificial galaxy. 
These spectra reproduce the characteristics of the NGC 
6166 data, assuming a range of surface brightnesses, 7j„ e 
{3.0, 1.0, 0.5, 0.1, 0.05} x I s ky, and sky subtraction errors, 
E in e {±0.02, ±0.01, ±0.005, ±0.001, 0.0} x I sky . The ar- 
tificial spectra were given a range of S/N ratios equivalent 
to summing several CCD rows as in §2.2.3 (10 < N rows < 
1000, depending on the desired S/N ratios). 

With the many variations described in the previous 
paragraph, 3 x 10 4 artificial spectra were created with 
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5 < S/N < 200 per A. The key results of fitting these sim- 
ulated galaxy spectra with a perfect template are shown 
in Figure 9. In the upper panels, for several ranges of 
"observed" velocity dispersions (a b s ), we show the frac- 
tional error in a, {(Jobs — fin/cobs, as a function of S/N 
ratio. If the data points cluster about the solid lines, then 
our methods are free from any systematic errors (template 
mismatch is not explored in this figure). The dashed lines 
show the best-fit curve of (<J bs—Vin/vobs oc [log (5 / N)]~ 2 . 
For (j b s < 700 km/s, the measured velocity dispersions are 
free from systematic errors for S/N > 10 per A. Figure 9 
suggests the presence of systematic errors when cr Q&s > 700 
km/s and S/N < 30 per A, but this is an artifact of the re- 
striction ai n < 900 km/s (in other words, a bs — 800 km/s 
with an error of —20% does not appear because there were 
no artificial spectra with <7j„ = 960 km/s). Based on the 
upper panels we conclude that the velocity dispersion pro- 
file derived in §4 is free from systematic errors (other than 
template mismatch, discussed below). 

The lower panels of Figure 9 show the formal errors, S bs, 
in a ba as a function of S/N ratio. The solid lines show the 
best-fit functions of the form 5 bs/&obs oc [log(SyiV)]~ 2 and 
illustrate the correlation of the estimated formal errors in 
a with S/N ratio. These estimates for the formal errors 
can be empirically verified using the scatter in the upper 
diagrams. The best-fit functions (using the same form), 
which represent the correlations of the true random er- 
rors with S/N ratio, are shown by the dashed lines in the 
lower panels. The agreement between the two curves is ex- 
cellent and indicates that we have satisfactorily estimated 
the random errors in a . 

The simulations were also used to explore whether fit- 
ting the Ca H & K features improves one's ability to re- 
solve large values of the Doppler broadening. The 4000A 
break and the H & K lines are broad features with presum- 
ably sufficient power on the scales required to accurately 
measure large velocity dispersions {a > 500 km/s). How- 
ever, incorporating that region of the spectrum does not 
improve the measurement of a because the S/N ratios are 
much lower than those in the region around the G band, a 
more useful feature because of the confluence of the higher 
fluxes and the intrinsic 2000 km/s width of the feature. 

We also used the simulations to explore the effects of 
template mismatch. By fitting the artificial galaxy spectra 
with spectra of individual late-type stars, we found small 
systematic errors on the order of < 5%. These systematic 
errors are not correlated with S/N ratio, and if one adopts 
the template that produces the lowest x 2 , the error is min- 
imized. Such systematic uncertainties, at a level of a few 
percent, are likely to be present in our velocity dispersion 
profile for NGC 6166. 

The simulations shown in Figure 9 also point to an addi- 
tional effect that will impact future studies: for a given de- 
sired fractional error in a, large velocity dispersions require 
data with higher S/N ratios. For example, measuring an 
internal velocity dispersion of a = 900 km/s with an accu- 
racy of ±10% requires spectra with S/N = 50 per A. An 
accuracy of ±15% requires S/N ~ 30 per A. Fortunately 
we did not find such large velocity dispersions in §4, even at 
a projected radius of 60 kpc. In the future, however, when 

8 The integrals are evaluated numerically using the ODEPACK and 
(http : //www . python . org/) scripts. 



the observer pushes towards larger radii and/or lower sur- 
face brightnesses in other targets, such levels of Doppler 
broadening may occur, and acquiring data with such high 
S/N ratios will be necessary to produce credible measure- 
ments of the internal kinematics. This requirement may 
be difficult to meet with current 8m-10m telescopes (even 
using the sky-subtraction technique of (Scmbach & Tonry 
1996). 

In summary we conclude that the velocity dispersion 
profile derived in §4 is free from systematic errors at the 
level of ±3%. Most of the systematic error arises from a 
mismatch of the template star with the galaxy. 

6. SIMPLE SPHERICALLY SYMMETRIC MODELS 

Under dynamical equilibrium, the observed velocity dis- 
persion profile is a constraint on the density profile (Jeans 
1915). In principle one can invert the observed kinemat- 
ics to directly infer the underlying distribution of mass 
(Binney & Mamon 1982) or use orbital-superposition 
approaches to build realistic mass models (Schwarzschild 
1979; Richstone & Tremaine 1984; Rix et al. 1997). How- 
ever, the construction of anisotropic mass models requires 
measurements of the Gauss- Hcrmite moments ft. 3 and /14, 
and we did not achieve sufficiently high S/N ratios to mea- 
sure these higher order moments of the LOSVD in the 
outer parts of NGC 6166 ( > 40 per A is required; van der 
Marel & Franx 1993), and our measurements in the inner 
parts were corrupted by the interloper (§4). 

With kinematics from only one position angle, and with- 
out the high-order moments of the LOSVD, several as- 
sumptions must be made before deriving meaningful con- 
straints on the mass profile. For the modeling of NGC 
6166 and Abell 2199, we assume spherical symmetry and 
zero net rotation. Even though the isophotes of NGC 6166 
have ellipticities ranging from 25% at R = 5 kpc to to 
50% at R — 75 kpc (Fig. 8, and the cD halo appears 
to be rotationally flattened, this rotation accounts for less 
than < 5% of the kinetic energy in the cD halo. Because 
neither the bulk of the cD, nor the cluster member kine- 
matics, is supported by rotation, the mild rotation in the 
outskirts of the cD has little impact on the modeling. We 
also note that the cllipticity of the isophotes does not nec- 
essarily invalidate our assumption of spherical symmetry 
- the total gravitational potential is most likely close to 
spherical given the morphology of the cluster. 

Thus this section is devoted to fitting simple spherically 
symmetric 1- and 2-component mass models to the veloc- 
ity dispersion profile of NGC 6166. The reduction of the 
Jeans equations to spherically symmetry greatly simpli- 
fies the computation of the projected line-of-sight velocity 
dispersion profiles, <j(R), for the model mass profiles we 
explore below. 

For non-rotating spherical systems the Jean's equations 
(Equations 4-55, 4-57a, and 4-60 of (Binney & Tremaine 
1987) are: 8 

^^ = -^(^ + ^ + 2/?(r)) (5) 
r y d\nr dlnr ' 

r°° vr A r 

QUADPACK (http://www.netlib.org/) routines from within Python 
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Fig. 10. — The observed velocity dispersion profile of NGC 6166 and its halo from integrated starlight (filled circles) and Abcll 2199 from 
cluster member kinematics (crosses). The solid line is the projected line-of-sight velocity dispersion profile of the best-fit isotropic single- 
component stellar mass model in which a power-law stellar density profile is used. The dashed line illustrates that the rising velocity dispersion 
profile of NGC 6166 can be well fit by an anisotropic model, while the right-hand panel shows the anisotropy used in the two models. Such 
high levels of anisotropy are not seen in other giant ellipticals and BCGs (typically —0.5 i> j3 ii 0.5; (Saglia, Kronawittcr, Gerhard, & Bender 
2000; Kronawitter et al. 2000; Gerhard ct al. 2001), and result in mass models that cannot support the velocity dispersion of the cluster 
unless the galaxies also follow such strongly tangcntially anisotropic orbits. 
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where (3{r) = 1 — Vg/v 2 , v is the stellar density at r as con- 
strained by the surface photometry, and M(r) is the total 
gravitational mass enclosed within radius r. The gravita- 
tional mass, M(r), includes the stellar mass and relies on 
assumptions about the radial gradient in the stellar M/L 
ratio. With the radial gradients in M/L B suggested in §3, 
the stellar density decreases more quickly than in the case 
of constant M/Lb- However, in the case of a decreasing 
M/Lb gradient, the light-weigted integrals in Equations 
6 and 7 are weighted towards larger radii, where the stars 
are increasingly more luminous per unit mass. 

For each parameterized model, we perform a search for 
the parameters that minimize the square of the residuals 
between the model and observed velocity dispersion pro- 
files. We choose to uniformly weight the data because 
weighting by the inverse of the formal errors would cause 
the minimization algorithm to find the parameters that 
best fit the inner regions of the galaxy. Furthermore, we 
exclude the inner R = 2" = 1.1 kpc from the fit because: 
(1) this region is heavily affected by the seeing (recall: 
0"8 FWHM); and (2) the spectrum in the center shows 
evidence for nuclear activity (also noted by Fisher, Illing- 
worth, & Franx 1995) and it is likely that a massive black 
hole dominates the kinematics in the very center, given the 
correlation between galaxy and central black hole masses 
(Gebhardt et al. 2000). 

Given the simplicity of the models, we attempt no 
derivation of uncertainties in the parameters of each model 
as their meaning could only be as valid as any given model 
is an appropriate representation of the mass distribution, 
or as valid as any of the underlying assumptions. Instead 
we use the \ 2 minima to assess whether given mass mod- 
els are consistent with the data at a given confidence level 
(> 95%). 

6.1. Single- Component Mass Models 

In this section we fit the data with models containing a 
single mass component. Two single-component mass mod- 



els are considered here: (1) a distribution of stars alone; 
and (2) a distribution of dark matter alone, in which the 
stars are effectively massless tracers of the dark potential 
(the stellar M/L is identically zero). While fitting single- 
component mass models is unlikely to provide physically 
meaningful characterizations of the distribution of mass 
in the core of Abell 2199 given the presence of both stars 
and dark matter, these models provide natural starting 
points for illustrating to what extent additional compo- 
nents might be required. After a brief discussion of the de- 
ficiencies of the single-component models, two-component 
mass models will be discussed in §6.2 

6.1.1. Single Component Stellar Model 

For this exercise, the observed velocity dispersion profile 
of NGC 6166 is fit using a mass density model defined by 
the fit of a King (1966) model to the surface brightness 
profile (sec §3, under the assumption of a constant M/L* 
ratio. The structural parameters for this model are listed 
in Table 1, leaving the overall normalization of M/L* and 
the anisotropy as the only free parameters. For this model 
we explore two variants: the isotropic case and the case in 
which the anisotropy, (3(r), has been adjusted to provide 
a reasonable match to the dispersion profile. 

We show these models in Figure 10(a), using the solid 
line to denote the best-fit isotropic model and the dashed 
line to illustrate the best-fit anisotropic model. Figure 
10(b) shows the anisotropy required to reproduce the ob- 
served velocity dispersion profile. These models illustrate 
that the stellar mass alone cannot support the rising a(R) 
profile in the cD halo, unless the stellar orbits become 
significantly (and unreasonably) tangcntially anisotropic 
(also see Tonry 1983). Such levels of anisotropy are in- 
consistent with what has been observed at 1-2 r e in other 
brightest cluster galaxies (typically —0.5 < [3 < 0.5; Saglia, 
Kronawitter, Gerhard, & Bender 2000; Kronawitter et al. 
2000; Gerhard et al. 2001). 

The purely stellar models imply a total mass of M* = 
4 x 10 12 M Q . This mass is insufficient to support the ob- 
served velocity dispersion of the cluster galaxies (TA2199 = 
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Fig. 11. — The observed velocity dispersion of NGC 6166 and its halo from integrated starlight (filled circles) and Abell 2199 from cluster 
member kinematics (crosses) as a function of projected radius with the single-component dark matter models overlayed. The solid, dashed, 
and dotted lines show the best-fit a = 0, 1, 1.5 gencralizcd-NFW density profiles, respectively. 



775 ±50 kms , unless the galaxy kinematics become even 
more strongly tangentially anisotropic (essentially follow- 
ing an extrapolation of the dashed line in the figure). Fig- 
ure 10(a) shows the velocity dispersion profile of the clus- 
ter derived from the 127 member galaxies within R < 1.5 
Mpc (obtained from NED). Given that the orbits of clus- 
ter galaxies are typically close to isotropic (Natarajan & 
Kneib 1996; van der Marel et al. 2000), it is unlikely 
that the galaxy orbits in Abell 2199 are as anisotropic 
as would be required with the purely stellar mass model. 
Even if M/L* is not a constant, under the constraints 
of §3, M(R) converges even more rapidly. As a result, 
the outskirts of NGC 6166 would be required to be more 
strongly anisotropic than in the case of a constant M/L*. 

We consider the single-component stellar models unfea- 
sible because both the stellar and galactic orbits are re- 
quired to be unreasonably tangentially anisotropic in order 
to support the large observed velocity dispersions. 

6.1.2. Single Component Dark Matter Model 

We now model the data using a single-component dark 
mass profile, in which the stars of NGC 6166 are massless 
tracers of the potential. In order to compute the projected 
line-of-sight velocity dispersion profile for NGC 6166, we 
use the stellar luminosity density profile (the projected 
velocity dispersion for the galaxy is a %/it-weighted mo- 
ment of the LOSVD). At larger radii, where the kinemat- 
ics of cluster galaxies are sampled, we assume that the 
galaxy number density traces the density of the dark mat- 
ter model. This assumption is valid because van der Marel 
et al. (2000) found that galaxy number density profiles of 
the CNOC1 clusters traced the mean dark matter density 
profile well. 

In Figure 11 we show the best-fit models using 
generalized- NFW density profiles with a6 {0,1,1.5}. Be- 
cause the NGC 6166 data-points have smaller relative un- 
certainties, the cD data carry more weight than the cluster 
data, and the extremely cuspy a = 1.5 model simply can- 
not be adjusted to fit both the inner regions of the cD while 
maintaining the overall normalization of the total mass as 
constrained by the cluster data. Even for the softer ha- 
los, a single component model cannot fit the galaxy and 
the cluster velocity dispersion profile simultaneously. The 
soft a < 1 model fits the rise in a(R) and the cluster data 



quite well, but leaves a hole in the inner region where the 
non-zero stellar mass of the cD resides. 

The concentration parameter, c = i?20o/ r s, is left as a 
free parameter in the fit. The best- fit a = 1 profile has 
c = 4 and i? 2 oo = 1-6 Mpc. Bullock et al. (2001) discussed 
the scatter in concentration parameters as a function of 
halo mass, and while they do not address halos as massive 
as ours (M 2 oo = 5 x 10 14 M Q ), an extrapolation of their 
Figure 4 indicates that c = 4 would fall within the ±lcr 
range of concentration parameters at this halo mass. 

Figure 11 also illustrates an important physical effect 
in the mass modeling: that if the density profile of the 
cD halo does not trace the number density profile of the 
galaxies in the cluster, then the velocity dispersion profiles 
of the two components will not be identical at a given pro- 
jected radius, (also see Natarajan & Kneib 1996). This 
effect appears as a discontinuity between the model veloc- 
ity dispersion profiles of the cD and the cluster. As the 
minimization algorithm adjusts the shape of the dark mat- 
ter density profile, the inferred number density profile of 
cluster galaxies may not have the same shape as the stel- 
lar density profile in the outskirts of the cD. As a result, 
the best-fit dark matter model may imply a distribution of 
galaxies that traces the potential differently than the stars 
in cD halo do (at a fixed projected radius). This differ- 
ence results in a velocity dispersion profile for the cD that 
may not seamlessly join the velocity dispersion profile of 
the cluster. Under different assumptions (forcing the stars 
in the cD halo to have originated from the current distri- 
bution of cluster galaxies), the velocity dispersion profiles 
of the cD and the cluster can be made a single, continu- 
ous function. Perhaps a future analysis of the color and 
absorption line strength gradients will help to constrain 
nature of the relationship between the stars in the cD halo 
and the cluster galaxies. 

In this section we restricted ourselves to isotropic dark 
matter models, a reasonable assumption because the 
CNOC1 clusters appear to be isotropic (van der Marel 
et al. 2000), and also because X-ray observations of Abell 
2199 constrain the anisotropy in the dark matter velocity 
distributions to —0.5 < (3 < 0.5 over a wide range of clus- 
ter radii (Mahdavi 2001). The conclusions drawn from 
the models in this section are insensitive to these levels of 
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Fig. 12. — The observed velocity dispersion of NGC 6166 and its halo from integrated starlight (filled circles) and Abell 2199 from cluster 
member kinematics (crosses) as a function of projected radius, with simple two-component models overlayed in which (a) a power-law density 
profile or (b) King (1966) model is used to represent the stellar mass. These models assume a constant B-band stellar M/L* ratio and the 
solid, dashed, and dotted lines indicate the models with the gcncralizcd-NFW dark matter halos for which a = 0,1,1.5, respectively. The 
models with a = have stellar M/L ratios M/Lb = 9, consistent with the M/L* ratios seen in other giant ellipticals (Gerhard et al. 2001). 
With a = 1, M/Lg = 1.2 and 0.2 when the two stellar mass density profiles are used, respectively. Using a = 1.5 forces M/Lg = 0, and 
these models still have far too much dark matter in the core to adequately match the observed velocity dispersion profile. 



anisotropy. 

6.2. Two- Component Mass Models: Stars + Dark Matter 

The previous section showed that one-component 
isotropic mass models do not satisfactorily reproduce the 
projected line-of-sight velocity dispersion profile of NGC 
6166. With this discrepancy in mind, we take the dark 
models and add a component that physically represents 
the stellar mass of the cD. The stellar mass profile is di- 
rectly constrained by the surface photometry, with the 
overall normalization (the stellar M/L ratio) to be con- 
strained by <j(R) in the inner regions of the galaxy. 

In this section we find the best-fit two-component 
isotropic mass models (see earlier references to Saglia, Kro- 
nawitter, Gerhard, & Bender 2000; Kronawitter et al. 
2000; van der Marel et al. 2000; Mahdavi 2001, regard- 
ing isotropy). The two mass components are expected to 
represent the stellar mass of the cD and the dark matter 
distribution of the cluster. Because the surface brightness 
profiles of giant ellipticals are well-fit by power-law and 
King (1966) density models (King 1978; Lauer 1985; Gra- 
ham et al. 1996), we explore models in which the stellar 
component is modeled cither by a power-law mass density 
profile or by a King model density profile, with structural 
parameters given in Table 1. The dark matter distribution 
of the cluster is parameterized by the generalized-NFW 
halo profile. 

In Figure 12(a), we show the best-fit models for the 
case where a power-law density profile of constant M/L* 
ratio is used to parameterize the stellar mass of the cD 
and where generalized-NFW density profiles with a e 
{0, 1, 1.5} are used to parameterize the dark matter halo. 
Figure 12(b) uses a King (1966) model for the stellar mass 
profile, and the dark matter halo density profiles are pa- 
rameterized as in 12(a). When fitting models with an 
a = 1 halo, we find concentration parameters of c = 3.5 
and 4.5 depending on whether the stellar mass is modeled 
as a power-law or King model, respectively. As above, 
these c values are consistent with the ±la range for the 
concentration given by Bullock et al. (2001). 



Of these two-component models, the dark matter halo 
profiles with a = fit best, yielding a B-band M/L* ra- 
tio of M/L b = 9 for either parameterization of the stellar 
density profile. If a = 1, the power-law density and King 
paramctcrizations have M/ Lb = 1-2 and 0.2, respectively. 
Using a dark matter halo with a = 1.5 cannot reproduce 
the observed velocity dispersion profile because the halo 
density profile has too much dark matter in the core (thus 
forcing M/L* = 0). The stellar M/Lb ratio derived using 
a = is consistent with what has been derived for other gi- 
ant ellipticals and BCGs (8 < M/L B < 12 for M > lO n M ; 
Gerhard et al. 2001). The a — 1 model is unlikely because 
the deduced M/Lb is inconsistent with the B — R color of 
the galaxy (for realistic initial mass functions, an extrapo- 
lation of the (Vazdekis et al. 1996) models to M/Lb = 1 
and B — R= 1.5 implies a mean stellar population age of 
24 Myr and [Fe/H] = +2.4). 

Allowing for the gradient in M/Lb, as constrained in 
§3, does not solve the problem. As stated earlier, the color 
gradient implies that M(R) converges rapidly for the stel- 
lar component. While the implication is that there is less 
stellar mass at large radii, the normalization of the stellar 
mass density profile is determined by the velocity disper- 
sions within R < 10 kpc. Together these effects produce a 
defneiency of stars at large radii, but there is no net im- 
pact on the model velocity dispersion profiles because the 
dark matter halo is adjusted during the % 2 minimization 
to maintain the total mass profile of the system. Further- 
more, if the stellar mass profile converges very rapidly, as 
is the case where the color gradient arises from a gradient 
in the mean stellar population ages, more dark matter is 
required to support the motions of the stars in the out- 
skirts of the galaxy. This need for additional dark matter 
at such radii is found by shrinking the core radius of the 
dark matter halo, but not by making the core cuspier (see 
earlier discussion). As a consequence, our conclusions are 
actually insensitive to the assumption of a constant M/L* 
ratio. Only in systems without a massive, compensating 
component of dark matter at large radii, will the gradi- 
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ent in stellar populations have a significant impact on the 
inferred total mass profile. 

6.3. Additional Points 

There are some additional limits to our data that should 
be noted. For example, the two-component mass models 
that employ the power-law density profile for the stellar 
mass produce a steeper a(R) gradient beyond R = 20 
kpc than the equivalent models that use King stellar mass 
profiles. This difference arises because the power-law den- 
sity profile has infinite extent, with the density falling as 
p oc r 2 - 5 ~ 3 , and the King model has finite extent. The 
surface photometry is not well-fit by an isothermal sphere 
and, as a result, the King model fit is truncated at a finite 
radius. In the integration to large radii there will sim- 
ply be more light at larger radii with a power-law stellar 
density profile and thus a steeper gradient in <r(R). 

We also investigated whether a King model was a bet- 
ter parameterization for dark matter density profile. Both 
the King model and a = generalized-NFW halo have 
similarly shaped density profiles in the inner regions, but 
they transition, respectively, to p oc r -2 and p oc r~ 3 at 
large radii. Unfortunately the current data do not extend 
beyond 1.5 Mpc, and we are therefore unable to determine 
whether the density profile of the dark matter halo declines 
either as p oc r -2 or as p oc r -3 . A better determination 
of the cluster density profile at large radii awaits member- 
ship surveys that extend beyond i?2oo (e-g-, as m Rincs 
et al. 2000). Such data would not only provide greater 
leverage on the shape of the mass profile at large cluster 
radii, but such membership information would allow for 
a test of the assumption that the galaxy number density 
traces the dark matter density. 

6.4. Comparison with X-Ray Masses 

At large radii, the inferred mass of the system is insen- 
sitive to the model adopted for the stellar density profile 
and mildly sensitive to the adopted dark matter density 
profile. With the a = generalized-NFW halo, we obtain 
M{< 0.5ft." 1 Mpc) = 1.7 x 10 14 /i- 1 )M Q (ft, = H /100) 
when adopting either the power-law or King model stellar 
density profile (and M{< 0.5ft," 1 Mpc) = 1.5 x 1O 14 /i _1 M 
with a = 1). These mass estimates are consistent with 
the mass derived from ASCA and ROSAT X-ray data: 
M(< 0.5ft" 1 Mpc) = (1.45 ±0.15) x 10 u h- 1 M e (Marke- 
vitch ct al. 1999). 

For radii within ~ 0.1ft, -1 Mpc, most X-ray mass de- 
terminations have been hampered by the relatively poor 
spatial resolution of the instruments and by what may 
be a central cooling flow (Siddiqui, Stewart, & John- 
stone 1998; Markevitch et al. 1999). As a result, there 
has been some variation in the published mass estimates: 
Markevitch et al. found M(< 0.1ft -1 Mpc) = (0.33 ± 
0.06) x 10 14 ft -1 M o , whereas Siddiqui et al. (1998), who use 
only the ROSAT PSPC data, found M(< 0.1ft -1 Mpc) w 
0.22 x lO 14 ft -1 Af (from their Figure 15). The difference 
between Markevitch et al. and Siddiqui et al. stems from 
the need to extrapolate the mass profile into the region of 
the cooling flow: Markevitch et al. (1999) modeled the X- 
ray temperature profile of the cluster using a polytrope 
and Siddiqui et al. (1998) modeled the cluster using a 
King profile. Because polytropes are, by definition, cus- 



pier than the isothermal cores of King models, Siddiqui 
et al. (1998) infer a shallower mass profile. Using X-ray 
data from Chandra, with much higher spatial and spectral 
resolution, Johnstone et al. (2002) suggest that the NFW 
profile provides a marginally better description of the total 
gravitational mass profile than a = 1.5 or a non-singular 
isothermal sphere (also see §6.1.2). Using their best-fit 
NFW model, those authors infer M(< 0. lft -1 Mpc) w 
0.15 x lO 14 ft _1 M0. Despite the uncertainties in the X- 
ray mass determinations, their estimates are in reasonable 
agreement with the values yielded by our two-component 
models: M(< 0.1ft" 1 Mpc) = 0.17 - 0.19 x 10 14 ft -1 M Q . 

6.5. Reconciling Soft Halo Profiles with Strong 
Gravitational Lensing 

Given its low redshift, there are currently no gravita- 
tional lensing data for Abell 2199. However, it is possible 
to ask whether the strong lensing characteristics of our 
our two-component mass models are consistent with the 
properties of known strong lenses at intermediate redshifts. 
Figure 13 shows a plot of the tangential critical radii of our 
two-component mass models when placed at the redshifts 
ziens S {0.2,0.3,0.4,0.5}. We restrict the figure to only 
include the two-component models with dark matter halo 
density profiles with a = and 1, because these models are 
able to reproduce the internal kinematics of NGC 6166 and 
the surrounding cluster. We mark the radii for tangential 
arcs in several intermediate redshift clusters, using crosses 
when the source redshift is known, and horizontal dotted 
lines when the source redshift is unknown. The horizontal 
shaded bar marks the i?-band effective radius (see Table 
1 that NGC 6166 would have at zi ens . At zi ens = 0.2, the 
strong lensing would be difficult to observe because the 
resulting arcs would be deep inside the image of the cD. 
Beyond zi ens > 0.3, the tangential critical radii lie well out- 
side r e and any strong lensing would be observable. The 
radii of arcs in clusters at zi ens > 0.3 are in rough concor- 
dance with our model predictions. 

As an additional note, the model with an a = NFW 
halo yields larger critical radii than the a = 1 NFW halo. 
This effect occurs because a more massive central galaxy 
is needed to reproduce the observed velocity dispersions 
within the central 10 kpc, and this extra structure in the 
mass profile gives the lensing signal a boost. While the 
current modeling has neglected departures from spherical 
symmetry and substructure, such effects serve to increase 
the predicted critical radii and improve the agreement 
(Bartelmann et al. 1995). Furthermore, if the dark matter 
halo density profile transitions at large radii to p ~ r -2 
instead of p ~ r" 3 , the excess surface mass densities can 
increase the predicted critical radii by factors of a few. 

In Figure 13, the predicted critical radii increase with 
zi e n S for a given z source . However, there appears to be no 
significant correlation of the observed arc radii with zi ens . 
The natural distribution of 9 cr it may not be as narrow 
as indicated by the observed arc radii, but may be trun- 
cated at low 6 cr it by the presence of large brightest cluster 
galaxies and by the limited spatial resolution of ground 
based imaging. These effects are likely to be the cause of 
the poor agreement between the predicted critical radii of 
our mass models at zi ens — 0.2 and the arc radii seen in 
z s=a 0.2 clusters. 
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Fig. 13. — The projected radii for the tangential critical curves as a function of source redshift for several possible lens redshifts for our 
two-component mass models. We show the best-fit two-component mass models in which the stellar density profile is parameterized by a 
power-law and the dark matter halo is parameterized either by a = (thick dashed) or by a = 1 (thick solid). The thin dashed and thin solid 
lines show the critical radii derived from the dark matter halos of those two-component models alone, ignoring the mass of the cD. Using 
the compilation of tangential arcs in Williams et al. (1999), we show example arcs seen in intermediate redshift clusters, using crosses when 
the source redshift is known and horizontal dotted lines when the source redshift is unknown. The horizontal gray bars indicate the effective 
radius NGC 6166 would have at each zi ena . If Abell 2199 were redshifted to cosmological distances, the cluster's mass profile would produce 
arcs similar to those seen in intermediate redshift clusters. For the case of the a = dark matter halo, the stellar mass serves as a large boost 
to the lensing signal. In the two-component mass model with a = 1, the stellar mass-to-light ratio is too low for the addition of the stellar 
mass to significantly affect the critical radii. Note that at z; ens = 0.2, the surface brightness of the cD would likely mask any strong arcs. 
Furthermore, because the greatest probability for lensing occurs when z source = 2z; ens , the predicted 6 cr it would also be difficult to observe 
with the limited spatial resolution of ground-based imaging. 



The strong gravitational lensing by MS2137-2353 has 
been used to constrain the mass profile of the cluster 
(Miralda-Escude 1995), and to predict a rising line-of- 
sight velocity dispersion profile for the central galaxy. The 
predicted velocity dispersion profile for the central galaxy 
of MS2137-2353 is similar to the velocity dispersion profile 
of NGC 6166, and in modeling the lensing cluster, Miralda- 
Escude (1995) also ruled out dark matter density profiles 
with p a r _1 in the inner region of the cluster. 

To summarize the mass profile we derived reproduces 
our observation of the kinematics in NGC 6166 and the 
surrounding cluster and is also consistent with the X-ray 
mass determinations for Abell 2199 and with observations 
of strong lensing in intermediate redshift clusters. 



7. DISCUSSION 

The cD galaxy NGC 6166 clearly shows a rise in a with 
projected radius, nearly reaching the observed velocity dis- 
persion of the galaxies in the surrounding cluster. NGC 
6166 is only the second cD galaxy with a measured veloc- 
ity dispersion profile that rises at large radii, after IC 1101 
(Dressier 1979). A rising velocity dispersion profile is not 
seen in the simulation by Dubinski (1998), in which he con- 
structs a cD-like galaxy through the continuous accretion 
of rotationally-supported spiral-like galaxies. 

The soft (a < 1) uniform density core of NGC 6166's 
dark halo is also inconsistent with cosmological iV-body 
simulations with inert CDM (Dave et al. 2001). Many 
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authors have now begun to explore other forms of cold 
dark matter, such as self-interacting, warm, or differently 
mutated forms of dark matter (Spergcl & Steinhardt 2000; 
Yoshida et al. 2000; Dave et al. 2001; Colin ct al. 2000; 
Dalcanton & Hogan 2001). The most extreme cases of 
active dark matter produce flat (a = 0) halo cores, but 
some fine-tuning of the dark matter interaction cross sec- 
tion may be required to avoid creating a universe of per- 
fectly spherical halos (Miralda-Escude 2002). Perhaps 
the discrepancy will be resolved with a being a function 
of scale (Jing & Suto 2000); though see (Klypin et al. 
2001). 

For the moment, the simulations are difficult to inter- 
pret. Not only is the scatter in a poorly known, but 
the quoted values of a have been derived by fitting the 
halo profiles outside r > 0.03i?2oo 9 (Klypin et al. 2001), 
and the consequences of the implied extrapolation of the 
parameterized models to r — > are ignored. Taylor & 
Navarro (2001) argue that the phase-space density profiles 
of dark matter halos are well-approximated by power-laws 
over 0.01i?2oo S r < -R200, and those authors extrapolate the 
power-law phase-space density profile to r = to obtain 
a = 0.75 (for NGC 6166 a = 0.75 implies M/L B = 4.5, 
too low by a factor of two, Gerhard et al. 2001). In- 
cluding the effects of baryons will likely be important on 
these scales as the cooling of baryonic matter will result in 
adiabatic contraction (Blumcnthal, Fabcr, Flores, & Pri- 
mack 1986; Rix et al. 1997; Keeton 2001), which steepens 
the inner profile of the halo (Alvarez, Shapiro, & Martel 
2001). It has also been suggested that the transfer of an- 
gular momentum from the baryons to the dark matter can 
soften the cores of halos (Weinberg & Katz 2001; Klypin, 
Zhao, & Somerville 2002). Furthermore, Weinberg (2001) 
suggests a number of dynamical processes that can drive 
the cores of dark matter halos towards a w 0. 

Future simulations will hopefully resolve these problems 
either (1) by using more physically motivated boundary 
conditions in the extrapolation to r — > (cf. Taylor & 
Navarro 2001), or (2) by resolving the physics within 
r < 3%i?2oo- Because Weinberg (2001) and Nusser & Sheth 
(1999) suggest that the inner profiles of dark matter halos 
may be sensitive to their initial conditions and histories 
of collapse and merging, the simulations should ultimately 
(1) predict the cosmic scatter in a (at a given halo mass), 
and (2) provide the means to estimate the timescales of 
formation for dark matter halos of a given mass and a (in- 
stead of using the concentration parameter; Bullock et al. 
2001). 

8. CONCLUSIONS 

We use the stellar kinematics of NGC 6166 and its halo 
to probe the mass distribution in the core of the rich clus- 
ter Abell 2199. The projected velocity dispersion along 
the major-axis of NGC 6166 rises steeply with radius be- 
yond the projected radius R w 10 kpc (where /ib ~ 23 
mag/arcsec 2 ). The velocity dispersion profile decreases 
from about 300 kms -1 at R = kpc to 200 kms -1 at 
R = 2 kpc (recall, the seeing was 0"8 FWHM = 0.45 
kpc), remains flat to R w 10 kpc, and then steadily rises 
to 660 ± 120 kms _1 at R w 60 kpc (where fj, B ~ 26 
mag/arcsec 2 ), nearly reaching the velocity dispersion of 

9 In our models 0.03i?200 ~ 50 kpc! 



the cluster (CTA2199 = 775 ± 50 kms -1 ). There is also evi- 
dence for mild rotation in intraclustcr stars beyond R — 20 
kpc ( I V/ cr I = 0.3), which is not seen in the kinematics of 
the cluster members. Because we measure the velocity dis- 
persion gradient with high precision, we have constrained 
the mass profiles of both the stellar and dark matter distri- 
butions, under the assumptions of isotropy and spherical 
symmetry. 

The primary results of our work indicate that the grav- 
itational potential of the cD halo is dominated by dark 
matter outside R — 20 kpc. Models in which the dark 
halo is represented by an a = 1.5 generalized-NFW density 
profile (e.g., Moore et al. 1999) are not able to reproduce 
the observed dispersion profile. While the a = 1 NFW 
profile is an improvement, the implied M/L B — 1 is not 
consistent with what has been measured for other giant el- 
liptical galaxies (8 < M/L B £ 12; Gerhard ct al. 2001). In 
the models that utilize the a = generalized-NFW dark 
matter halo, one obtains M/Lb = 9, a more reasonable 
value. 

The two-component mass models are consistent with the 
mass profiles derived from X-ray observations (Markevitch 
et al. 1999; Siddiqui, Stewart, & Johnstone 1998) over 
the range of radii for which our observations overlap. We 
also derived tangential critical radii for our two-component 
mass models assuming Abell 2199 were placed at cosmo- 
logical rcdshifts. The predicted tangential critical radii are 
consistent with the radii of arcs in strong lensing clusters 
at intermediate redshift. 

Our results obtained for Abell 2199 have strong impli- 
cations for the study of dark matter halo profiles. The 
best-fit models are inconsistent with the predictions of N- 
body simulations of cold dark matter. However, we can- 
not yet verify the universality of dark matter halo profiles. 
Simulations arc needed to provide more robust estimates 
of a at radii r < 0.01r2oo (see, e.g., Klypin et al. 2001) 
and to predict the scatter in a. Observationally, we are 
obtaining kinematical information for additional cD galax- 
ies with the goal of providing empirical constraints on the 
scatter in dark matter halo properties. 

With the possibility of observing cD galaxies to large 
radii using very long 10') slits, the task of background 
subtraction will become easier and more accurate. Us- 
ing a large database of cD galaxy kinematics, with sev- 
eral long-slit position angles per galaxy, we can explore 
more complicated models, in which both intrinsic flatten- 
ing and anisotropy are included. At the same time, the 
expansion of catalogs of cluster galaxies (e.g., Rines et al. 
2000), and dramatic increases in the spatial resolution of 
X-ray data, will enable us to better constrain the inner 
and outer slopes of dark matter halo profiles, while allow- 
ing us to avoid many of the simplifying assumptions. Such 
improvements may soon lead to the confirmation or refu- 
tation of universal dark matter density profiles, at least 
for the largest dynamically bound systems. 
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Table 1 

Structural Parameters for the Stellar Profile 11 







r 1 / 4 -law 6 




King Model 






Power-Law d 


Filter Correction 6 r e (kpc) 


rms (mag) 


r c (kpc) 


ri/2 (kpc) 




rms (mag) 


r c (kpc) 


7 


rms (mag) 


B 


None 


61.3 


0.073 


3.84 


102 


10.25 


0.071 


2.35 


1.22 


0.034 


R 


None 


47.5 


0.058 


3.58 


46.2 


9.50 


0.069 


2.45 


1.26 


0.048 


B 


Age 


19.7 


0.064 


2.82 


18.2 


8.22 


0.033 


2.94 


1.50 


0.059 


R 


Age 


19.7 


0.052 


2.78 


18.2 


8.23 


0.029 


2.92 


1.49 


0.070 


B 


Fe/H 


34.2 


0.064 


3.30 


31.6 


8.85 


0.042 


2.66 


1.35 


0.046 


R 


Fe/H 


34.2 


0.053 


3.28 


31.6 


8.86 


0.051 


2.63 


1.35 


0.055 



Note. - ■ a To avoid the effects of the seeing disk, we used an inner fitting radius of rj„ = 2"5. We also restricted the 
outer fitting radius to r out — 75 kpc because that radius is also the outer radius for the stellar kinematics. b The r 1 / 4 -law 
has been fit so that the resultant half-light (effective) radii can be compared to the half-light radii derived using the other 
profiles. c The King (1966) profile has two free parameters: r c , the core radius, and ^o/a 2 , the truncation energy (see Binney 
& Tremaine 1987). d The power-law, of the form p oc [1 + (r/rc) 2 ] 1 , has two free parameters as well. Note that because the 
best-fit power-law profile has 7 < 2, the total mass/luminosity is divergent and there can be no half-light radius. e The color 
gradient may be interpreted as a spatial variation in the stellar population ages and/or metallicities (see text), and we use 
this interpretation to explore corrections to the radial variation in (M/L*). Even though the color gradient in NGC 6166 is 
small (see text), the inferred half-mass radii are very sensitive to the correction. 
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Table 2 



Major-Axis Kinematics for NGC 6166 



R 


AR 


a 


5a 


V 


5v 


x 2 


(kpc) 


(kpc) 


(kms- 1 ) 


(kms -1 ) 


(kms -1 ) 


(kms -1 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


-17.16 


6.29 


404.7 


47.5 


-40.4 


34.1 


1.36 


-12.47 


2.61 


326.1 


59.1 


14.0 


27.1 


1.21 


-9.74 


2.37 


267.1 


26.6 


-12.3 


15.8 


1.14 


-7.66 


1.31 


243.6 


23.7 


-12.7 


13.7 


1.07 


-6.35 


0.83 


207.1 


26.6 


-19.8 


15.1 


1.17 


-4.69 


2.02 


219.9 


13.5 


-10.1 


7.6 


1.49 


-2.97 


0.95 


216.7 


12.4 


-10.9 


7.0 


1.30 


-1.96 


0.59 


200.3 


14.5 


0.9 


7.5 


1.45 


-1.25 


0.36 


215.4 


14.4 


4.2 


7.9 


1.27 


-0.71 


0.24 


264.9 


14.3 


0.7 


8.9 


1.18 


-0.24 


0.24 


305.7 


18.8 


-2.6 


5.5 


1.35 


0.24 


0.24 


289.0 


16.9 


0.9 


6.7 


1.27 


0.71 


0.24 


247.2 


15.3 


-3.9 


9.7 


1.44 


1.25 


0.36 


209.5 


13.1 


4.5 


8.3 


1.39 


1.90 


0.47 


208.2 


13.7 


-9.2 


8.2 


1.29 


2.73 


0.71 


217.9 


14.2 


4.1 


7.8 


1.18 


3.80 


0.24 


226.4 


15.0 


-18.7 


7.9 


1.04 


3.86 


1.07 


216.2 


15.4 


6.7 


8.4 


1.20 


4.27 


0.47 


240.8 


10.1 


-7.1 


5.3 


1.77 


4.93 


0.59 


220.1 


9.3 


-2.8 


4.6 


1.46 




1 ^/L 
1 .0^ 




10. D 


1.0 


7 8 


1 1 fi 

1 . 10 


5.58 


0.47 


197.3 


12.0 


4.0 


5.0 


1.34 


6.17 


0.47 


191.6 


11.8 


17.1 


4.8 


1.22 


7.01 


0.95 


209.5 


11.3 


4.0 


4.2 


1.48 


8.01 


3.21 


233.9 


19.6 


6.1 


7.6 


1.12 


8.19 


1.19 


213.0 


11.8 


-5.6 


4.2 


1.49 


9.44 


1.07 


221.7 


14.8 


-10.4 


5.5 


1.36 


10.69 


1.19 


282.7 


15.6 


8.1 


6.6 


1.41 


12.05 


1.31 


262.9 


16.1 


-27.7 


8.0 


1.18 


13.42 


1.42 


248.0 


16.8 


5.3 


10.6 


1.36 


15.67 


2.85 


275.0 


15.3 


-4.7 


9.4 


1.37 


18.64 


2.85 


316.1 


17.4 


-5.3 


11.7 


1.25 


27.49 


2.97 


351.1 


29.3 


38.1 


19.5 


1.39 


30.81 


3.44 


445.6 


60.5 


45.0 


37.9 


1.35 


34.91 


4.75 


490.9 


58.6 


53.1 


29.9 


1.72 


39.30 


4.04 


529.1 


73.9 


114.3 


47.3 


1.32 


57.77 


24.34 


663.0 


121.2 


232.6 


61.8 


2.38 
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